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Metallic nanostructures are of immense scientific interest owing to unexpectedly strong interaction with 
light in deep subwavelength scales. Resonant excitations of surface and cavity plasmonic modes mediate 
strong light localization in nanoscale objects. Nevertheless, the role of surface plasmon-polaritons (SPP) in 
light transmission through a simple one- dimensional system with metallic nanoslits has been the subject of 
longstanding debates. Here, we propose a unified theory that consistently explains the controversial effects 
of SPPs in metallic nanoslit arrays. We show that the SPPs excited on the entrance and exit interfaces induce 
near-total internal reflection and abrupt phase change of the slit-guided mode. These fundamental effects 
quantitatively describe positive and negative effects of SPP excitation in a self- consistent manner. 
Importantly, the theory shows excellent agreement with rigorous numerical calculations while providing 
profound physical insight into the properties of nanoplasmonic systems. 

Renewed by the discovery of extraordinary optical transmission (EOT) 1,2 , extensive study has been devoted 
to explaining light transmission through metal films with subwavelength aperture arrays. Initially, EOT 
through hole arrays was understood by local field enhancement with interfacial excitations evanescently 
coupled through subwavelength holes. Surface plasmon-polaritons (SPPs) in the optical domain 3 " 5 and geomet- 
rical surface resonances in the THz or microwave spectral ranges 6 induce such coupled interfacial excitations. 

However, subwavelength slit arrays have shown many distinguished behaviors from those of hole arrays. A 
widely accepted enhancement mechanism for slit arrays is based on Fabry-Perot resonances of slit-guided 
modes 7 " 18 , i.e., cavity modes (CMs). Many recent papers conclude that CMs provide the enhancement mechanism 
while SPPs play only a negative role 13 " 16 . In this view that deters SPPs, spectral location of the transmission 
minimum corresponds to the SPP resonance condition and the associated field pattern shows typical SPP 
character with its null at the aperture opening. Strong surface-plasmonic absorption 15,16 , excitation of a non- 
resonant SPP 13 , and the surface-plasmonic bandgap effect 1718 have been suggested explanations for these negative 
effects of SPPs. The negative effect of an SPP has also been reported for hole arrays, and different viewpoints based 
on the nonresonant SPP excitation 19 and destructive interference between transmission pathways via surface 
modes and hole-guided modes 20 have been published as suggested explanations. 

Nevertheless, another bundle of recent papers have reported enhanced transmission with clear SPP characters 
in surface field patterns and frequency- dispersive properties 7 " 12 . In these analyses, SPPs enhance transmission 
with coexisting CMs 10 " 12 . Associated with SPP-CM hybrids, transmission peaks become much narrower than the 
pure CM resonances 8 " 10,12 ' 1718 — the CM resonance condition shifts abruptly 8,12 and an asymmetric Fano profile 
appears 21 " 23 . Therefore the role of SPPs in metallic nano aperture arrays is still under debate, and the essential 
physics remain unclear due to these diverging interpretations. 

In this paper, we show that controversial SPP-related effects can be consistently described by a single unified 
model that treats a metallic nanoslit array as an optical cavity with SPP-resonant boundaries. We theoretically 
prove that all aforementioned SPP-related effects such as the antiresonant extinction, null field at the aperture 
opening, bandwidth narrowing, and abrupt shift of the CM resonance condition are rooted in a single resonance 
interaction: a surface-plasmonic Fano resonance that occurs when the external light and CM couple at the 
entrance and exit interfaces. Contributed by the SPP excitation, metal-film interfaces act as a Fano-resonant 
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Figure 1 | Schematic of a surface-plasmonic Fano resonance model of a 
metallic nanoslit array, (a), Schematic of a one-dimensional array of 
metallic nanoslits. (b), Coupling of external light with the fundamental 
CM at the entrance and exit interfaces, t represents the coupling coefficient 
between the external light and CM while p in denotes the internal reflection 
coefficient of the CM. (c), Fano-resonance interpretation of the coupling 
between the external light and the CM at the single interface. p SP and t sp 
represent reflection and transmission coefficients via an SPP state, 
respectively, while p D and t d denote the nonresonant reflection and 
transmission amplitudes, respectively. rj in (rj ex ) is the coupling probability 
of the SPP with the CM (external light). 

gate that closes or opens nanoslit cavities and causes associated phase 
changes in the internal reflection of the CM. This interfacial inter- 
action successfully describes various metamorphic SPP -related 
effects in a physically intuitive manner. 

Analytic theory 

Consider transmission of transverse-magnetic (TM) polarized light 
through a metal film perforated by an array of slits with period A, 
thickness d, and slit width w as shown in Fig. la. For deep sub wave- 
length slits (w <C X) in an optically thick metal film (d > skin depth 
y27T|£ M '| 1/2 )> the light transmission can be described by a Fabry- 
Perot formula for the fundamental CM 7 ' 8,12 ' 24,25 as 



T = 



l- Pin V^ 



where ft = ft' + is the complex propagation constant of the CM, t 
is the coupling coefficient between the CM and external planewave, 
and p in is the internal reflection coefficient of the CM as illustrated in 
Fig. lb. Equation (1) is generally applicable to cases with arbitrary 
angles of incidence as long as only zero -order waves propagate in the 
surrounding media. Enhanced transmission peaks appear when the 
multiple scattering denominator in Eq. (1) becomes minimal at the 
phase-matching condition P'd q = (q + l)7i-arg(p in ), where q is an 
integer. The role of an SPP is implicitly held in t and p in as a mech- 
anism causing a Fano resonance to occur at the top and bottom 
interfaces. 

Here, we further examine the coupling processes at each film 
interface. We treat each interface as a Fano-resonant boundary 
where an SPP acts as a discrete state that interferes with the non- 
resonant continuum. In Fig. lc, an SPP-resonant pathway interferes 
with a nonresonant pathway in the scattering processes of the incid- 
ent CM. An SPP originally excited by the incident CM emits the 
reflected CM with probability f] in and the transmitted external 



radiation with probability f] ex . This process provides the resonant 
components of the single-interface reflection p SP and transmission 
t sp coefficients. The incident CM at each interface also couples to the 
non-plasmonic reflection p D and transmission t d . Using the optical 
Fano resonance theory developed by [22] , the single-interface trans- 
mission t and reflection p in coefficients are written as 



t(<5) = td + 



1-iS 

= ?SP 



Pm($) = Pv + 




(2) 



(3) 



where 3 = (co — <x> SP )/2y tot is the normalized frequency, co SP is the 
resonance frequency of the SPP, y tot is the total decay rate of the SPP, 
?] in and ?] ex represent radiation probabilities of the SPP to the CM and 
the external planewave, respectively, (j) is the plasmonic transmission 
phase at 3 = 0, and (p = arg(p SP ) is the plasmonic reflection phase at 3 
= 0. In this view to treat the interface responses, the time reversibility 
requirement dictates the phase differences between the plasmonic 
and non-plasmonic contributions, i.e., £ = arg(i D ) — arg(T SP ) and 
£ = arg(p D ) — arg(p SP ) at 3 = 0, to be determined by 22 



e it = _ ^radl T D| 

IpdI 



+ i\ 1- 



^rad 2 !^! 2 



(4) 



(5) 



where rj md = rj in + r] ex is the total radiation probability of the SPP. 
Once establishing the real- valued parameters such as radiation prob- 
abilities (rj in and rj ex ) and magnitudes of the nonresonant pathways 
(\pr>\ and |t d |), we can describe the interference between the coup- 
ling pathways and the spectral properties of t and p in in a fully 
deterministic manner with these phase relations. 

Lossless and nondispersive system 

We check the theory with numerical calculation based on the rig- 
orous coupled- wave analysis 26 (RCWA) and finite element method 27 
(FEM). The metal is modeled by a complex dielectric constant £ M = 
£ M ' + *W> where s M '< 0 and £ M "^ |£m'|- Figures 2a and 2b show 
film- thickness -dependent transmission spectrum T(k,d) due to the 
RCWA and our analytic theory in Eqs (1) ~ (5), respectively. We 
assume here a lossless (e M " = 0) metallic slit array with £ M = — 5 and 
slit width w = 0.05A as an ideal case that reveals the essential physics 
with minimized complexity. We will discuss more realistic cases of 
lossy dispersive metals later in this paper. We extract the SPP res- 
onance parameters cl> sp , y, f/ in , and f] ex for the analytic theory cal- 
culation from the single-interface SPP excitation spectrum calculated 
by FEM (see method section for details). We note that no numerical 
fitting method was used to find these parameters. Figures 2a and 2b 
confirm excellent quantitative agreement of periodically appearing 
CM Fabry-Perot resonance peaks; more notably, the figures confirm 
drastic modification of CM resonance properties involving peak 
extinction at XI A =1.13 (red dashes) and an abrupt peak shift over 
the bright background region around XI A = 1.06 (black dashes). 
These CM resonance modification effects are deeply associated with 
responses of the film interfaces and are therefore central to the con- 
troversy over the role of SPPs. 

We analyze the response of the interface by assuming a semi- 
infinitely-thick slit array. Figure 2c shows the average intensity of 
the surface-normal electric field \E Z \ 2 at the metal-air interface under 
planewave incidence from air (red squares) and CM incidence from 
slit (blue squares). The average \E Z \ 2 spectra at the interface exhibit 
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Figure 2 | Film transmittance T(a and b) vs. optical response of the single interface (c ~ f). Thickness-wavelength (^-d) map of transmittance (T) 
through a finite film due to (a), RCWA and (b), the analytic theory. Note wavelength X and film thickness d are in the unit of period A. (c), SPP excitation 
spectra (\E Z \ 2 at the interface) for light incidence from air (red, □) and guided-mode incidence from slit (blue, □) on the single interface, (d), \t\ 2 
(black), I t sp 1 2 (blue), and I t d 1 2 (red) spectra. Open square symbols ( □ ) are numerically calculated by FEM while curves ( — ) are obtained by the analytic 
theory, (e), Phase difference between t sp and t d due to Eq. (4). (f), Internal reflection phases arg(p in ) due to FEM ( □) and Eq. (5) ( — ). e M = —5, w = 
0.05A, and surface-normal incidence are assumed for all cases. 



symmetric Lorentzian profiles with a common resonance center X SP 
= 1.062 A and full- width at half-maximum (FWHM) A^ S p = 
0.03211 A. We attribute this surface excitation to a pure SPP on a 
patterned surface. Using a pure-SPP model developed by Liu and 
Lalanne 5 , we predict the SPP resonance wavelength and bandwidth 
(FWHM) on a metallic slit array to be 

^SPF 



As- 



A^sp = As - 



1- arg(fs + r s )/27i' 



^sp , , _ i x 2nsp /r 



SPF 



-ln(>s + rs| 



"SP 



(6) 



where t s and r s are in-plane SPP transmission and reflection coeffi- 
cients at a single isolated slit, respectively, ^ SPF = n SP 'A is the SPP 
resonance wavelength on a flat metal surface, and n SP + n SP " = [s M / 
(1 + £m)] 1/2 is the complex effective index of the SPP on a flat metal 
surface. FEM calculation of SPP scattering by a single isolated slit 
yields arg(f s + r s ) = -0.1047ti and |f s + r s | _1 = 1.1097. The SPP 
resonance wavelength and bandwidth due to Eqs. (6) and (7) are X SP 
= 1.0624A and A^ SP = 0.03346A. These values quantitatively agree 
with those obtained from the surface excitation spectrum in Fig. 2c. 
See Supplementary Section I for the derivation of Eqs. (6) and (7) and 
Supplementary Section II for the FEM calculation of the in-plane SPP 
transmission t$ and reflection r s coefficients. 

The pure SPP excitation and associated Fano resonance at the 
interface successfully describe the drastic modification of CM res- 
onance properties. The single-interface transmittance \z\ 2 due to the 
analytic theory (solid curve — ) in Fig. 2d shows excellent agreement 
with the numerical calculation result (square symbols □) due to 
FEM. A typical Fano resonance profile appears with its resonant 
enhancement peak at X R = 1.058 A (<^sp) an d antiresonant extinc- 
tion at X AR = 1.1 30 A (>A, SP ). First, the destructive interference 
between the plasmonic (t sp ) and non-plasmonic (t d ) coupling coef- 
ficients explains the antiresonant extinction of CM resonance peaks 
at X AR . Figure 2e shows the spectral behavior of the phase difference 
between t sp and t d , £ = arg(T D )-arg(T SP ) due to Eq. (4). Note at X = 
X AR the phase difference £ = n where |t S p| — |t d | as shown in Fig. 2d, 
leading to complete destructive interference followed by z = 0. The 
null field at the slit opening for the excitation by external planewave 
in Fig. 3a is a natural consequence of the null excitation of the 
CM due to the complete destructive interference at X AR . In this 
interference description, the null field at the slit opening does not 
require any special electromagnetic excitation such as nonresonant 



SPP 1319 . Note that in Fig. 3b the field at the slit opening is non-zero 
for the interface excitation by a CM at X = X AR . The effect of the 
antiresonant extinction of z on the film- transmittance T is obvious: 



the internal reflectance of the CM becomes total ( 



iPii 



1 



1) as t approaches 0. Therefore, the slits at X = X AR behave as closed 
cavities that generally support vanishingly narrow and extremely 
high-quality cavity resonance peaks. 

On the other hand, at the SPP enhancement condition for X = X SP 
or X R , \z\ 2 is maximal and consequently the slits behave as open 
cavities with partial transmission and reflection at the interface as 
shown in Figs. 3c and 3d. We note strong SPP excitation at the 
interface that leads to resonant phase change in the internal 
reflection of the CM. In Fig. 2f, the internal reflection phase arg(p in ) 
shows an S-shaped 2n change centered at X = X SP . This phase beha- 
vior clearly describes the abrupt shift of the CM resonance peaks 
observed in Figs. 2a and 2b. Recalling the CM resonance condition 
P'd q = (q + l)n — arg(p in ), 2n phase-change in p in results in a 
transition of the resonant film thickness d q —> d q + 2 . This transition 
in the CM resonance condition is widely found in the literature 811 ' 28 
but has not been explained in terms of internal reflection phase 
change associated with the surface-plasmonic Fano resonance. 

Lossy system 

Now we consider the effect of material dissipation. Figures 4a and 4b 
show transmission spectra for several different film thicknesses (d/A 
= 1.03, 1.06, and 1.1) in lossless (s M = — 5) and lossy (e M = — 5 + 
0.01/) cases, respectively. All solid curves calculated by the analytic 
theory in Eqs. (1) ~ (5) are again in excellent agreement with the 
RCWA simulation results (square symbols □ ). For the lossy case in 
Fig. 4b, the material dissipation is included in the analytic theory by 
including the non-radiative decay rate of the SPP y nr = (2nd 
^sp)("sp"/n S p') ( see Supplementary Section I for details) and the 
complex propagation constant of the CM ft = fi' + ifi" that is given 
by the equation 29 g M tanh[w(^ 2 - k 0 2 ) m l2] = - [(/J 2 - 8 M k 0 2 )/(P 2 - 
/c 0 2 )] 172 , where k 0 = 2n/X and c is speed of light in vacuum. 

The transmission peaks for q = 4 in Figs. 4a and 4b clearly reveal 
the effect of material absorption on the CM resonance peaks near the 
antiresonant extinction condition X = X AR = 1.1 30 A. For the lossless 
case in Fig. 4a, the linewidth of the transmission peak tends to vanish 
as the peak approaches ^ A r- Diverging CM localization lifetime in 
the closed cavity regime leads to an extremely narrow linewidth and a 
high local field enhancement. This is pointed out as the origin of the 
diverging inter-slit coupling matrix at the surface-mode resonance 
condition in the modal expansion method 2 ' 9 . A transmission peak 
with an extremely narrow linewidth quickly disappears with material 
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Figure 3 | Magnetic field (H y ) distributions for single interface excitation for (a), planewave incidence from air at the antiresonance condition 
(X = ^ar)> and (b), guided-mode (CM) incidence from slit at X = X AR , (c), planewave incidence from air at the SPP resonance center (X = X SP )> 
and (d), guided-mode (CM) incidence from slit at X = X SF . 



dissipation while those with a relatively wide linewidth near X = X SP 
= 1.062 A or X R = 1.058 A are loss-insensitive as shown in Figs. 4a 
and 4b. 

The relation between the peak transmittance (T max ) and the res- 
onance quality factor Q, i.e., the number of effective oscillations in 
the cavity, quantitatively explains a narrow peak's high sensitivity to 
material dissipation. From Eq. (1) and the surface-plasmonic absorp- 



tion A SP = 1 



|Pm| 2 > we obtain 47i 2 T max ~ Q 2 |t| 2 and Q 



□ RCWA Analytic Theory 
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Figure 4 | Effect of material dissipation on the film transmittance ( T) and 
resonance quality factor. The film transmittance (T) through (a), lossless 
(e M " = 0) and (b), lossy (e M " = 0.01) metal films for d/A = 1.03 (blue), 
1.06 (black), and 1.1 (red). Square symbols (□) are obtained by RCWA 
simulation, and curves ( — ) represent the analytic theory based on Eqs. ( 1 ) 
~ (5). (c), Resonance Q-factor spectra for e M " = 0 (blue) and 0.01 (red), 
(d), Ratio of the SPP-induced absorption (A SP ) to the CM-induced 
absorption (A CM ) in a lossy metal film (e M " = 0.01). ^ R = 1.058 A and 
X AR = 1.1 30 A. 



27iA SP " 1 exp( — P"d q ) in a high Q-factor regime (Q ^> 1). Whereas 
T max = 1 in the lossless case, it decreases with Q 2 as ohmic damping 
of free electrons causes surface-plasmonic absorption (A SP ) at the 
interface and propagation loss of the CM (A CM ) inside the cavity. In 
Fig. 4c, the decrease in Q for the lossy case (s M " = 0.01) is negligible 
near X = X SP or X R where Q is relatively small. In the extremely high 
Q band near ^ AR , however, Q is remarkably suppressed for the lossy 
case (red curve), resulting in strong suppression of T max indicated by 
the gray curve enveloping the peaks in Fig. 4b. It is also worth noting 
that the surface-plasmonic absorption is not a dominant absorption 
channel responsible for the CM resonance extinction at ^ AR . The 
ratio of A SP to A CM in Fig. 4d is 0. 1 at ^ A r> and thereby A CM is nearly 
10 times stronger than A SP . In contrast, at X = X SP = 1.062 A where 
surface-plasmonic absorption is maximized (A SP /A CM = 7.44), the 
effect of loss on Q and T max is not remarkable as shown in Figs. 4b 
and 4c. Therefore, the propagation loss of a CM is dominantly 
responsible for the resonance extinction in this model system with 
s M = -5 + 0.01/. 

Lossy and dispersive system 

In our description, X R , X SP , and ^ AR are crucial parameters strongly 
dependent on the metal dielectric constant e M . We use FEM to cal- 
culate X R , X SP , and ^ AR as a function of the real part £ M ' of £ M , and the 
result is shown in Fig. 5. As £ M ' decreases to far negative values, X R 
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Figure 5 | Dependence of characteristic wavelengths of the surface- 
plasmonic Fano resonance on metal dielectric constant (<?m'): the 
antiresonant extinction (^ A r)> SPP resonance center (^ S p)> resonant 
enhancement (^ R ), and SPP resonance center on a flat, unpatterned 
surface (A, SPF ). We assume slit width w = 0.05 A. 
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Figure 6 | Film-transmittance T(a, d, and g), surface excitation \E Z \ 2 (b, 
e, and h), and single-interface transmittance | r| 2 (c, f, and i) spectra for 
silver slit arrays: a, b, and c for period A = 400 nm and slit width w = 
20 nm. d, e, and f for A = 500 nm and w = 25 nm. g, h, and i for A = 
400 nm and w = 40 nm. We use the silver dielectric constant (e M ) listed in 
[32]. Note that I tI 2 is on a linear scale, instead of the log scale used in 
Fig. 2d. 

and X AR approach the canonical Rayleigh anomaly ^R ay i e igh = A and 
the SPP resonance wavelength on a flat, unpatterned metal surface 
X SPF , respectively. Dependence of these characteristic locations on 
the metal dielectric constant (e M ') suggests that previous confusion 
in the role of SPPs may originate from the spectral proximity of X R 
and X AR to Vayieigh and ^ SPF . For example, at e M ' ~ — 15 where many 
previous analyses have been performed 2 , |^ AR — ^ SPF | ~ 4 X 10 " 4 A 
and |V — Vayieighl ~ 1.5 X 10" 3 A. In this situation, it is likely to 
form a hasty conclusion that the SPP plays only a negative ro le 15 ' 16 ' 30 ' 31 
and that the actual SPP-associated effects, such as the abrupt shift of 
the CM resonance condition, are confused with the effect of Rayleigh 
anomaly 28 . In [ 15,16] for example, the authors analyzed Au slit arrays 
in the near- and mid-infrared spectral domains where £ M ' < —30. 
Concluded by the coincidence of the antiresonant CM resonance 
extinction with the SPP resonance wavelength (^ S pf) on a flat, unpat- 
terned metal surface, they attribute the CM resonance extinction 
solely to the surface-plasmonic absorption. However, our model 
shows that an SPP induces high-Q CM resonances at the antireso- 
nance condition, and the CM resonances in this case are highly 
sensitive to losses in any kind including both surface-plasmonic 
and cavity-modal absorption. Moreover, it is evident in Fig. 5 that 
the true SPP resonance wavelength ^ SP on a perforated metal surface 
differs from both A, SPF and Vayieigh even down to £ M ' ~ —30, which 
corresponds to Ag for wavelength —800 nm. 



Finally, we show that all aspects of the simple Fano resonance 
model presented above are also present in realistic systems with lossy 
dispersive metals. In Fig. 6, we show the film-transmittance T(^,d), 
average field intensity |£ z | 2 , and single-interface transmittance \t(X) \ 2 
for silver nanoslit arrays with several different periods and slit 
widths. We use RCWA for T(X,d) and FEM for \E Z \ 2 and \t(X)\ 2 with 
the realistic 8 M of silver experimentally obtained by Johnson and 
Christy 32 . For all three different cases, we confirm the same char- 
acteristic features as for the simplified case in Fig. 2; they include CM 
resonance modifications at X = X SP and ^ AR in T(X,d) associated with 
a single Lorentzian surface excitation in the \E Z \ 2 spectrum and a 
typical Fano profile in |t(^)| 2 . In Figs. 6a ~ 6c for A = 400 nm 
and w = 20 nm, X AR is fairly separated from }t SPF and resonance 



shift (<L 



2 ) due to a 2n reflection-phase change that fully 



appears around ^ SP . In Figs. 6d ~ 6f for A = 500 nm and w = 
25 nm, no additional feature appears. Affected by larger |g M '| for 
the longer wavelength, X AR almost coincides with ^ SPF (\X AR — ^ SPF | 
~ 3 X 10" 3 A) and the resonance shift near X R is seemingly associated 
with the Rayleigh anomaly at X = A as previously discussed in Fig. 5. 
In Figs. 6g ~ 6i for A = 400 nm and a wider slit width w = 40 nm, 
the wider slit results in a wide bandwidth in the SPP excitation. The 
Rayleigh anomaly, defined as sharp intensity variations occurring 
when an evanescent higher-order wave turns into a propagating 
wave, causes a corresponding decrease in the zero-order intensity. 
Our model is unable to describe effects associated with the Rayleigh 
anomaly because it is limited to a subwavelength-period regime 
where only the zero-order waves are allowed in the radiation con- 
tinuum. Nevertheless, the rigorous calculation results for different 
cases in Figs. 1 and 6 consistently show no necessary effect of the 
Rayleigh anomaly on the resonance properties while the character- 
istic features of our interfacial Fano resonance model persistently 
appear. Therefore, the results in Figs. 1 and 6 suggest that the inter- 
facial Fano resonance is the fundamental origin of the cavity res- 
onance modification. In the previous literature, Sarrazin et al. 
reported a comprehensive spectral and surface field analysis also 
showing that the Rayleigh anomaly is unnecessary and they sug- 
gested the significance of surface-plasmonic Fano resonances with 
a phenomenological argument based on complex poles and zeros of 
scattering amplitudes 33 . 

We have shown that various aspects of SPPs in EOT through 
metallic nanoslit arrays can be consistently understood by the sur- 
face-plasmonic Fano resonance in the coupling of external radiation 
to the slit cavity mode. The Fano resonance interpretation was first 
suggested by Genet et al 21 in order to explain the asymmetric profile 
and red shift of the enhanced transmission feature. They assumed a 
single discrete state without any other localized states such as a slit- or 
hole-guided mode. With detailed coupling processes unclear, the 
original Fano resonance interpretation has been used for phe- 
nomenological analyses of experimental and numerical data 2 . Our 
theory clearly describes where the Fano -type interference occurs, 
how it modifies the optical response of a metal surface with periodic 
nanoslits, and how the SPP-resonant metal surface finally contri- 
butes to metamorphic cavity- resonance properties. 

Consistency with previous theories 

In addition, our theory provides deeper physical insight into the 
microscopic theory of EOT developed by Liu and Lalanne 5 . For 
two-dimensional hole arrays, they found the single-interface trans- 
mission coefficient 



t(co) = t d + 



2ab 



exp( - insvcoA/c) - (f s + r s ) 



(8) 



for normal incidence, where a and b denote coupling coefficients of 
an SPP with the hole-guided mode and external radiation, respect- 
ively. The Fano -type interference is an inevitable consequence of this 
expression as its first and second terms on the right-hand side 
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represent the nonresonant and resonant contributions, respectively. 
Equations (2) ~ (5) for Fano resonance in the single-interface coup- 
ling processes are applicable to two-dimensional hole arrays in prin- 
ciple. Therefore Eq. (8) should be consistent with our Fano resonance 
interpretation. Indeed, Eq. (8) reduces to Eq. (2) with radiation 
probabilities 

c \b\ 2 c \a\ 2 

^" |n S p|Ay tot |fs + rs| |n SP |Ay tot \t s + r s \ ' W 

and resonant transmission phase 

(j)= arg(a)+ arg(fc)- arg(f s + r s )- arg(« S p). (10) 

In Eq. (9), the total decay rate y tot = (nc/X S p 2 ) A^ SP with ^ SP and AX S p 
in Eqs. (6) and (7), respectively. See Supplementary Section III for 
detailed derivation. These relations describe how elementary scatter- 
ing processes of electromagnetic fields at the metal-film interfaces are 
associated with the more fundamental wave kinematic effect of Fano 
resonances. The formal consistency of our model with the micro- 
scopic theory of EOT suggests further importance of the Fano res- 
onance interaction in longer wavelength ranges beyond the visible 
domain. A series of theoretical 34,35 and experimental 36 analyses 
recently showed that an additional contribution from the quasi- 
cylindrical wave can also be described by the same elementary scat- 
tering coefficients a, b, f s , and r s of an SPP. In the near-infrared and 
longer wavelength domains, the quasi-cylindrical wave is known to 
significantly contribute to the resonances in periodic arrays of metal- 
lic nanoapertures 3637 . 

Conclusions 

In conclusion, we propose a surface-plasmonic Fano resonance the- 
ory of the light transmission through metallic nanoslit arrays. 
Importantly, seemingly paradoxical, metamorphic SPP-related 
effects are clearly explained by the pure surface-plasmonic Fano 
resonance effects at the film interfaces, which cause drastic modi- 
fication of the cavity- resonances inside nanoslits. We also show that 
the interfacial Fano resonance interpretation is formally consistent 
with the microscopic theory of EOT through two-dimensional hole 
arrays. Therefore, for a two-dimensional array of large holes that 
allow propagating guided modes, a surface mode must lead to fun- 
damentally the same effects on the cavity-mode resonance properties 
as those in one-dimensional slit arrays. For example, Catrysse and 
Fan 20 reported the antiresonant extinction of transmission peaks 
associated with hole-guided modes when the surface mode is reso- 
nantly excited in an SiC film with cylindrical holes. We believe that 
our theory unifies different interpretations and illuminates the origin 
of previous confusion regarding the role of SPPs. For example, the 
antiresonant extinction of CM resonance peaks is not simply a nega- 
tive SPP effect but is rooted in the SPP-induced total internal 
reflection of the CM (SPP-induced cavity closing). Note that, in this 
case, the SPP actually contributes in a positive way as it leads to very 
high-Q CM resonances. Therefore, appropriate loss compensation 
methods 38,39 are of great interest at the antiresonant extinction con- 
dition as extremely high-Q nanocavity resonances are expected. We 
summarize how our theory unifies previous partial interpretations in 
Supplementary Table 1. Our model is limited to the zero-order 
regime and deep subwavelength slits that allow only the fundamental 
guided mode. Further development of our approach to more general 
cases of interfacial Fano resonance coupling with higher-order prop- 
agating waves and multiple localized modes may yield deeper phys- 
ical insight into various nanophotonic and surface-plasmonic 
systems where interplay of coexisting modes induces versatile spec- 
tral properties and novel optical effects. 

Methods 

To estimate basic resonance parameters a> SP , y tot , rj in , and rj ex for the analytic theory, 
we use surface excitation spectra in Fig. 2c. Two excitation spectra (red squares) for 



planewave incidence from air (blue squares) and for CM incidence from slit are 
denoted by E ex and E in , respectively. E ex and E in show Lorentzian resonance peaks 
with a common center and bandwidth. First, eo SP and y tot are taken directly from the 
peak location and half- width at half-maximum of the Lorentzian profile. We obtain 
a» SP = 0.9419 X 2nc/A and y tot = 0.02847 X 2nd A (c is speed of light in vacuum). 
Second, rj in and n ex are taken from the peak values of E in and E ex . Considering Lorentz 
reciprocity theorem in the mode coupling processes, the radiation probability is 
proportional to the excitation probability. Therefore, rj in oc £ in (^ SP ) and rj ex oc 
^ex(^sp)- Including the relation for the total radiation probability n in + n ex = y ra d/y W> 
where y rad and y tot are radiation and total decay rate of the SPP mode, respectively, we 
obtain 

mn = £ in(^p) ?rad ^ ^ ^ex^Sp) y rad _ 

Em (^sp ) + E ex (Asp ) V tot (^sp ) + E ex (^sp ) V tot 

The expressions for y rad and y tot are given in Supplementary Section I. Using these 
relations, we obtain rj in = 0.1589 and rj ex = 0.8411 from Fig. 2c. In this calculation, 
identical incoming power is assumed for the two different cases of external planewave 
incidence and CM incidence. In addition, y rad = y tot in Fig. 2c as we assume lossless 
metal. 
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